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We explore the bulk-edge correspondence for topological insulators (superconductors) 
without time-reversal symmetry from the viewpoint of the index theorem for open 
spaces. We assume generic Hamiltonians not only with a linear dispersion but also with 
higher order derivatives arising from generic band structures. Using a generalized index 
theorem valid for such systems, we show the equivalence between the spectral flow of 
the edge states and the Chern numbers specifying the bulk systems. 
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^ ■ 1. Introduction 

q When a bulk gapped system has a topologically nontrivial ground state, gapless 

edge states appear if the system has boundaries. This bulk-edge correspondence has 

CN ! become increasingly important since the recent discovery of various types of topologi- 

> ■ 

cal insulators and superconductors EEIEEIE] For quantum Hall systems, the relationship 

between the Chern number of Landau levekP^ and the quantized Hall conductance as 

■ a result of the edge currentP was elucidated by HatsugaP for a tight-binding model. 

O ' 



Volovik^l and recently Essin and Gurarie^ formulated the correspondence by using the 
Green functions for more generic systems with a domain wall separating topologically 
distinct systems. An adiabatic approach was also developed by Sato et al. for topological 



superconductors.^ 

More generically, defects such as kinks, vortices, and monopoles can be regarded 
as a type of boundary, and indeed they allow midgap states. In particular, zero-energy 
states are known to reflect the topological property of the corresponding bulk systems. 
These zero-energy states were investigated by index theorems that show the equivalence 
between the analytical index concerning the zero-energy states as solutions of differential 
equations and the topological index associated with the topological configurations of 
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background fields included in the operators (Hamiltonians). Therefore, it is natural to 
utilize the index theorems to investigate the problem of edge states associated with the 
conventional boundary mentioned above. 

In this paper, we develop a method of investigating the bulk-edge correspondence 
for class A insulators and class D superconductors by using the index theorems formu- 
lated by Calliasp^ Weinberg,^ and Niemi and Semenoff J-^ and recently applied to the 
Majorana states of topological super conductors p^"^"^ which are valid for open infinite 
systems in arbitrary dimensions. Edge states are usually described as follows: Suppose 
that a translationally invariant boundary lies along the ?/-axis for simplicity. Then, we 
have a one- dimensional Hamiltonian with a parameter k y . Edge states are defined by 
the spectral flow of this Hamiltonian as a function of k y . To utilize the index theo- 
rems for such systems, one often regards the parameter of the Hamiltonian as a part 
of background fields. However, since the parameter in the present case, k y , does not 
have upper and lower bounds, it is difficult to calculate the topological index in a stan- 
dard manner such as by derivative expansions. Furthermore, the index theorems^'^E' 
studied thus far were restricted to the Dirac operator (Hamiltonian) only with the first- 
order derivatives, whereas in condensed matter physics, we often encounter systems 
with higher derivatives arising from general band structures. The purpose of this paper 
is to develop a generalized index theorem applicable to such generic Hamiltonians, and 
to present several techniques invented for its application to the bulk-edge problem. 

This paper is organized as follows: In §[21 we present the setup of the bulk-edge cor- 
respondence for a generic Hamiltonian. We first introduce a bulk Hamiltonian whose 
topological properties are controlled by one parameter. We next consider a Hamiltonian 
with a domain wall that separates topologically distinct systems, by making the parame- 
ter coordinate-dependent. This Hamiltonian is referred to as a domain wall Hamiltonian. 
For this Hamiltonian, edge states are expected at the domain wall. To characterize the 
edge states, we define the spectral flow of these states. In §21 we present the formulas for 
the Chern number characterizing the bulk systems. In §H we first extend the domain 
wall Hamiltonian, and show that the edge states can be regarded as the zero-energy 
states of this extended Hamiltonian. This implies that the edge states are characterized 
by the index of the extended Hamiltonian. In §51 we first derive a generalized index 
theorem and next apply it to the extended Hamiltonian. We show that the index is 
given by the Chern numbers, and, together with the fact that the spectral flow of the 
edge states is equal to the index of the extended Hamiltonian, we thus establish the 
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bulk-edge correspondence. 

2. Setup for Bulk-Edge Correspondence 

We investigate a generic bulk Hamiltonian 'H B (—id x , —id y ) depending on a constant 
parameter m controlling the topological properties of the bulk system. Without loss 
of generality, we assume that the gap is open at the zero energy and that negative 
energy states are occupied. We also assume translational invariance, which enables us 
to represent the Hamiltonian in the momentum space such that 

n B (k x ,k y ) = n B (k), (2.1) 

where k stands for the set of two momenta, k x and k y . Then, the occupied bands give 
rise to a nontrivial Chern (TKNN) number c(m), which shall be discussed in the next 
section. Correspondingly, we introduce a model with a domain wall at x = 0, which 
separates two different topological states, i.e., 

lim m(x) = m±. (2.2) 

x— »±oo 

Let us denote such a Hamiltonian as T-L(—id x , —id y , x), where the x dependence implies 
that the parameter m is a function of x. We assume that the Hamiltonian is still 
translationally invariant along the y direction in the presence of the domain wall, and 
hence, the Fourier transformation leads to the Hamiltonian 

H{-id x ,ky,x). (2.3) 

This is a one-parameter family of one-dimensional Hamiltonians and its spectrum as 
a function of the parameter k y includes edge states localized around the domain wall. 
The eigenvalue equation is 

H(-id x , k y , x)(f)j(x, k y ) = ej{k y )4>j{x, k y ). (2.4) 

Although the bulk system has a gap at the zero energy, eq. (12. 4p admits several states 
that cross the zero energy. These states are referred to as edge states. More precisely, 
the edge states are those with the spectrum £j(— Ajl)£j(Aju) < 0, where — A jL and Aju 
are the lower and upper bounds of k y within which £j(k y ) is defined as a continuous 
function.^ Assume that the model has N + edge states with Sj(Ajxj) > 0, and AL states 
with Ej(Aju) < 0. The spectral flow v is defined by 

v = N + - N_. (2.5) 
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The relationship between the Chern number c(m±) of the bulk system ( 12. ip and the 
spectral flow v of the domain wall model ( 12. 3 p is the generic bulk-edge correspondence 
investigated in this paper. 

In what follows, we suppose a Dirac Hamiltonian or a more generic Hamiltonian of 
a Dirac-like matrix structure with higher-derivative terms. For later convenience, we 
introduce the following notation separating the derivative terms from the others, 

H(-id x ,k y ,x) = ^ 7 i n) H^r + <^(^,x). (2.6) 

n>l 

We have assumed that m(x) is not coupled with the derivative terms. The coefficients 
of the derivative terms 7^ may be a function of k y . The k y dependence of does 
not change the following argument. 

A nontrivial example of the bulk Hamiltonian j^EIEIllIIlEIGSlllII %b _ _^ a .Q. _|_ 
cr 3 (m + bd 2 ), where j = 1,2, (xi = x,x 2 = y) and d 2 = d 2 + d 2 . Its momentum 
representation is 

U B (k) = a 3 k 3 + a 3 {m - bk 2 ), (2.7) 
where k 2 = k 2 + k 2 . The corresponding domain wall model is therefore 

H(-id x , k v , x) = -m x d x + a 2 k y + a 3 [m(x) - bk\ + bd 2 ] . (2.8) 
Comparing with eq. (I2.6p . we have 

l { x ] = a u lx ] = ~ba 3 , 

5U(k y , x) = a 2 k y + a 3 [m(x) - bk 2 y ] . (2.9) 
As shown in the next section, the bulk Hamiltonian (12.71) has a nontrivial Chern number. 

3. Chern Number of the Bulk System 

In this section, we summarize notational conventions for the Chern number char- 
acterizing the bulk systemP^m L e t, ip a (k) be the normalized eigenstate of T-L B (k) with 
the eigenvalue e a (k), 

n B (k)tfj a (k) = e Q (k)Mk), (3.1) 

where a specifies the band of the spectrum, and let ip(k) = fa, • • • ) be the set of 
wave functions with negative energies, s a (k) < 0. Then, 

P = ^ (3.2) 
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serves as a projection operator to the negative energy bands. The Berry connection 
one-form A and the curvature two-form F are defined by 

A = ^dtf), 

F = dA + A 2 = d^ t (l-P)d^, (3.3) 

where d is the exterior derivative with respect to the momentum fcj. The Chern number 
is defined by 

c = — / tr F 

= h e%1 S fd 2 k(d k M^l(d k ^ a ), (3.4) 

where in the second line, we have used the fact that 1 — P in eq. ( 13. 3 j) is the projection 
operator to the positive energy states. 
The following relation is also useful: 

P(dP) 2 = tpFtpK (3.5) 

It follows that trP = trP(dP) 2 . In the case of the model in eq. (12. 7p . this formula 
is useful, since the projection operator ( 13. 2 p can be expressed alternatively as P = 
\ (l — W B /R), where R(k) is defined by (H B ) 2 = R 2 t (R > 0) and given concretely by 
R{k) = [k 2 + (m — bk 2 ) 2 ] 1 ^ 2 . The Chern number is thus given by 



/ \ 1 f ,2i m + bk 2 

4vr J [jfea + ( m _ bk 2 ) 2 ] 3/2 

= g [sg n ( m ) + sgn(6)] . (3.6) 

Therefore, if m + and m_ in eq. (12.21) have opposite signs, a domain wall appears in model 
(12.81) that separates two distinct systems with trivial and nontrivial Chern numbers. 

4. Index of the Edge States 

Thus far, we have defined the spectral flow (12. 5p for the domain wall system as well 
as the Chern number (13.41) for the gapped ground state of the bulk system. To prove 
the correspondence between them, we introduce an extended Hamiltonian, 

H = -ivT X d ky + T 2 U(-id x , k y , x), (4.1) 
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where v (> 0) is a constant parameter. Although % is a one-parameter family of one- 
dimensional Hamiltonians, % is now a two-dimensional model. This technique was used 
in many references dealing with index theorems.^®^^ It should be noted that this 
Hamiltonian has a generalized chiral symmetry, {"H, t%} = 0. This enables us to regard 
the spectral flow of H as an index of the extended Hamiltonian ~H. 

To be specific, let us consider the zero energy states of this extended Hamiltonian, 

m {k y ,x) = 0. (4.2) 

Owing to the chiral symmetry of "H, these states can also be eigenstates of T3, and 
hence, we assume the following product wave functions: 

* ( I ) ' (4 ' 3) 

where <pj is the wave function in eq. (12.41) . and the former and latter in the above 
equation have r 3 = 1 and r 3 = —1, respectively. By applying an adiabatic deformation 
preserving topological invariants, we can choose <pj independently of k y . Then, it follows 
that 

=F ej) ^ = 0. (4.4) 

The solution of this equation is given by 

^.(^ oce^^^A (4.5) 

This tells us that for any positive value of v, {ipo ) is normalizable if and only if 
£j(— Aj L ) > > £j(Aju) (£j(— Ajl) < < ej(AjXj)). Therefore, we have 

-mdU = v. (4.6) 

The above argument clearly shows that the index is independent of v. This point plays 
a crucial role below. 

5. Index Theorem for the Extended Hamiltonian 

In the previous section, we introduced the extended model (14. ip whose zero energy 
states are given by the edge states of the domain wall system. We next apply the index 
theorem, which makes it possible to express the index as a topological invariant of the 
Hamiltonian. 
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5.1 Index theorem 

We will first generalize the index theorem for the present system including higher- 
order derivatives generically. Let us denote the extended Hamiltonian (14. 1|) simply as 

U = -iy.dj + 5U, (5.1) 

where j — 1, 2 and dj = d t] with (ti, t 2 ) = (fey, a;)- The 7 "matrices" and are defined 
by 



71 = OTi, 72 - 



n>l 

5U = T 2 5U(k y ,x). (5.2) 

The chiral symmetry of the Hamiltonian is denoted by {U, r 3 } = 0. The index of U in 
the previous section can be defined alternatively by 

mdU = lim Trr 3 — ■ — , (5.3) 

m->o iM-U 

where Tr stands for integration over t (= k y , x) as well as the trace in the r-space 

and the internal band structure. As has been shown in many references] El [T4l[T5 l[i6l[T7] 

above index can be expressed in terms of the divergence of a chiral current for the Dirac 

Hamiltonian. In what follows, we will show that even if the Hamiltonian includes higher 

derivatives, we can define a similar current. In the derivation of the index theorem, the 

formal identity^ 

iM U 

limtrr 3 ^5 2 (t - s) = limtv r 3 6 2 (t - s) + limtrr 3 ^5 2 (t - s), (5.4) 

s ^ iM - U ^ ^ iM-U 

plays a crucial role. We show in Appendix |~A] that the last term can be written as the 
divergence of a current 

U 1 

lin } trT3 ^ — T, 5 ^ ~ s ) = ~o 9 i J i> ( 5 - 5 ) 
iM -U 2 

where the current is defined by 

J At) = limtrr 3 r 7 - % -^S 2 (t - s), (5.6) 

iM-U 

with 

Ti = 71 = vn, 

r 2 = r^AH4-^). (5-7) 

n>l 
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Note that I^ includes the effect of the higher derivatives through a polynomial p n 
defined by 

p n (a,b)= Yl {-) l a n - l ~ l V. (5.8) 

0<i<n-l 

If the model has only the first-order derivative, p\(a,b) = 1 implies that T 2 becomes 
the usual matrix. The reason why we have introduced the polynomial (15.81) is to use 
the following identity: 

p n (a,b)(a + b) = a n - (-b) n . (5.9) 

This equation is required to derive eq. (15. 5p . as shown in Appendix I Al The index (15.31) 
can thus be rewritten as 

ind^ = lim TiT 3 f(U 2 /p 2 ) - - lim Tr&J,-, (5.10) 

/X— >oo 2 M— s>0 

where in the r.h.s., the first term is associated with the chiral anomaly when % is a 
Dirac Hamiltonian only with linear kinetic terms. To treat generic systems with or 
without higher derivative terms in a unified manner, we have introduced an ultraviolet 
regulator f(x), which is a smooth function satisfying /(0) = 1 and /(oo) = 0. 

In what follows, we shall calculate each term in eq. (I5.10p separately. In the con- 
ventional index theorem, the derivative expansion is applicable, since the background 
fields are assumed to be constant at infinities. In the present Hamiltonian (15. ip . the 
parameter m indeed satisfies this condition owing to eq. (12. 2p . However, 8% diverges in 
the limit \k y \ — > 00, as seen from eq. (12.91) . This implies that we cannot resort to this 
simple technique of derivative expansion. Nevertheless, a crucial role is played by the 
parameter v introduced in eq. (14. ip instead. As we have noted, the index (14 .6p holds for 
any value of v (> 0). This makes it possible to expand eq. (I5.10p with respect to this 
parameter. 



5.1.1 Calculation of the T3 term 

Let us first consider the T3 term in eq. (I5.10p . As a regulator function, we choose the 
standard function f(x) = e~ x . To calculate this term in the limit v — > 0, we explicitly 
write the integration over k y such that 

Trr 3 e-« 2 ^ 2 = / dk y [°° —e^Tr ^"^V^, (5.11) 
J J-oo 2tt 

where Tr^. stands for the trace Tr except for the integration over k y , and we have used 
the Fourier representation of the 5 function in eq. (15. 4p . 5(t 1 — si) = 5{k y — k' y ) = 
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/ ^e i( - ky and taken the limit k' y ->■ k y . Using 

and making the scale transformation u; — >• oj/v, we can obtain the following expansion 
with respect to v: 



-iuiky g— W 2 / giu)ky 



i--(T 3 d ky u-2iujd ky : 
a* 



-(to 2 +H 2 )/v 2 



0{v 2 



Thus, we obtain the expansion of the r 3 term with respect to v such that 



Tr r 3 e 



-n 2 /n 2 



dfc,, 



du 1 

27T f 



Tr x r 3 



1 ^{r 3 d ky U - 2iud ky 

t 1 



-^ 2 +H 2 )/^ 2 



(5.13) 



0{v 2 ). 
(5.14) 



In the above expansion, terms without r 3 vanish by the trace over r. Therefore, in the 
limit v — > 0, we have 



Trr 3 e 



-H 2 /n 2 



dh. 



°° duJ 1 Tr a .^c~ (cj2+K2)/ ^ 2 . 

.00 27T /i 2 ^ 



(5.15) 



where the factor of 2 originates from the trace over r, implying that Tr x here does not 
include this space any longer. The integration over u is carried out to give 



lim Trr 3 e 



-H 2 /» 2 



7T fi — ^00 



lim / dkyTr x 



i&u 

/i dk y 



(5.16) 



As we shall show, this term cancels out the next term below and does not contribute 
to the final index. 



5.1.2 Calculation of the Tr<9iJi term 

Let us first compute Ji(k y ,x) for any values of k y and x. To this end, we introduce 
the notation Q = (u,k x ) as conjugate "momenta" to t — (k y ,x). Then, by using the 
Fourier representation of the two-dimensional 5 function 5 2 (t — s) in (15.61) . we have 

J (2vr) 2 iM-H 

d 2 Q 1 , 

trOT 3 _ I , „. Q ZTnTTZ 7a Z H\ ' ( 5 - 17 ) 



(27r) 2 t x M + ivu + vd ky - T 3 H(k x - id x , k y , x 

where the minus sign in the last line originates from T\T 3 T\ = — r 3 . Rescaling u — y u/v 

and expanding with respect to v, we see that eq. (15.171) becomes 

f d 2 Q 1 
J 1 (k y ,x) = - / — —^trr 3 - — ■ — ^j + Oiy) 

J (27T) 2 IbJ + TiM - T 3 H 

f d 2 Q H , 

^ 2 J(W b ^TW' (5 ' 18) 
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where in the second line, the pref actor of 2 originates from the trace over the r space, 
and the limit M — > has been taken. At this stage, we arrive at 

H 



lim / dx Ji(k v , x) = [ — Tr a 



= viK)- (5-19) 

Here, the function rj is a formal expression of the spectral asymmetry, and a certain 
suitable regularization is required. We adoptP^ 

2 f 00 

Vs(k y ) = -=j- ( — — prTTyT / dm/Tr^e"^ 2 . (5.20) 



r(( s + i)/2) 7 

In the limit s — > 0, ?7 S — > 7] = rj as calculated above It then follows thatP^ 

r duu*d u (Tr x u^Le-^ 2 ). (5.21) 



dfc y r((s + l)/2) J \ dfc y 

Since this should be well-defined, we can safely take the limit s — > 0, and hence, we 
obtain 



limTr^^e-"^ 2 . (5.22) 



y/n u->o dk 
Thus, the Tr<9i</i term in (15.101) can be integrated as 

i r „ , /■ 9 



2 



y d 2 t9iJi = y d^y^T y dxJi(A; 2/ ,x 



dfc 



d7 7( fc y) 



lim / d^Tr^-^e"^^ 2 , (5.23) 

i^OO J " /idfcy 



where we have set u = Therefore, we see that the leading term of Tr<9iJi with 
respect to v is the continuous part of the spectral asymmetry of "H. It turns out that 
this term (15.231) cancels out the r 3 term (15.161) . 



5.1.3 Calculation of the Trd 2 J 2 term 

Because of the cancellation of the two terms above, the remaining Trd 2 J 2 term is 
solely involved in the index. Carrying out the integration of the d 2 J 2 term over t 2 = x 
yields 

jdHd 2 J 2 = j dky[J 2 (ky, OO)- J 2 (ky, -OO)]. (5.24) 
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Therefore, we need to compute the current J 2 itself at the spatial infinities. In the 
limit of t 2 = x — > ±00, we have assumed m{x) — > m± = const. Then, we can replace 
H(—id x , k y , x) by H B (-id x , k y ), 

H -)• -ivr x d ky + r 2 U B (-id x , k y ). (5.25) 

Then, we have 

e^ nt Ue mt = -in [ivu + vd ky - T 3 H B (k x - id x , k y )] 

= -in [ivu + vd ky - r 3 n B (k x , k y )] + 0(d x ). (5.26) 

Here, it should be noted that d x is irrelevant, since m is constant. Using these and 
Appendix [B[ we have 



J 2 {ky, ±00) 



d^ 

(27T) 



-trl T 



1 



dk x ) ivu + vd kv - r 3 H B 



(5.27) 



Now let us make the scale transformation u — >■ u/v and take the limit v — > 0. Then, 
the series expansion with respect to u yields 



^(/c^, ±00) 



d 2 Q 



27T 



,-tr 1 T 



1 



1 



1 



d 2 l] 



rtr 



V iu — T^H B iu — T^H B V itU — T-yhL B 

iu 1 „ 1 



+ 0(v) 



V U* 



B\2 



-o, 



IU 



U B y iu-H 



B 



+ 0(v), (5.28) 



where the trace in the r space has been taken in the second line. The first term, being an 
odd function of u, vanishes, whereas the second term remains finite. It can be written 
using the Green function of the bulk system, 

1 



G(iu, k x , k y 



such that 



r d 2 n 

J 2 {k y ,±oo) = 2i J 7^p tr 



iu - H B (k x , k y ) 

m B ' 



(5.29) 



d. 



2i 



d 2 Q 

W) 



rtr 



dk x ) iu-H B y iu-H 
dn B \ n dG 

Gr- 



dk 3 



dk. 



(5.30) 



As we will show in Appendix l~C| this becomes the Chern number c(m±) after the 
integration over k y . Thus, we reach the final expression 

/oo 
dky J 2 (k y , ±00) = 2c(m±). (5.31) 
-OO 
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5.2 Result for the index 

Thus far, we have calculated each term in the r.h.s. of the index theorem (I5.10p . 
Equations (15. 23|) and (15.311) allow us to write the index as 

1 f dr] 

indU = Trr 3 - - J — - [c(m+) - c(m_)j , (5.32) 

where the second term in the r.h.s. containing the integral represents the continuous 
part of the spectral asymmetry rj. This is analogous to the family index theorem derived 
by Niemi and Semenoff As shown in the calculations above, the first and second terms 
become eqs. (I5.16P and (I5.23p . respectively, and the index finally becomes 

— ind"H = c(m + ) — c(m_). (5.33) 

It then follows from eq. (14. 6 j) that the spectral flow of the edge states is equal to the 
difference in the Chern numbers, 

v = N + - iV_ = c(m+) - c(m_). (5.34) 

Thus, we establish the bulk-edge correspondence valid not only for the Dirac Hamilto- 
nian but also for a wider class of Hamiltonians with higher-order derivative terms. 

6. Summary and Discussions 

We have explored the bulk-edge correspondence for topological insulators and super- 
conductors with broken time reversal symmetry. We have shown that the spectral flow 
of the edge states, which are located near a domain wall separating two topologically 
distinct bulk states, is identical to the difference in the Chern numbers characterizing 
these two bulk states. To show this, we have utilized the index theorem developed by 
Callias, Weinberg, Niemi, and Semenoff, which is valid for open spaces. To this end, we 
have utilized several techniques. One is the dimensional extension of a one-dimensional 
Hamiltonian describing the edge states to a two-dimensional Hamiltonian, which was 
often used for index theorems in previous studies. Another is the introduction of the 
parameter v in eq. (14.11) in the dimensional extension. This allows us to expand the 
topological index with respect to this parameter. Since we have assumed Hamiltoni- 
ans including in general not only the linear kinetic terms but also those with higher 
derivatives, the proof of the correspondence in this paper is valid not only for the Dirac 
Hamiltonian but also for generic nonrelativistic Hamiltonians. 

The generalization of the present approach to the bulk-surface correspondence of 
three-dimensional topological insulators and superconductors, as well as to various topo- 
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logical classes other than Chern insulators will be interesting. As to a lattice model, 
Hatsugai has shown the bulk-edge correspondence for the tight-binding model on a 
square lattice in a uniform magnetic field by using the transfer matrix method. It will 
be interesting to develop a lattice version of the present method applicable to generic 
lattice models. 
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Appendix A: Detailed calculation of the divergence of the current 

In this Appendix, we calculate the r.h.s. of eq. (15. 5p . It follows from definition (15. 6|) 
that the divergence of the current is 



djJjit) = dUimtr t 3 T 



iM-n 



S 2 (t - s) 



limtrr^d' + d*) ^S 2 (t - s), 

s-+t ] J 'iM-n 



(A-l) 



where the repeated index j should be, of course, summed over. In what follows, we will 
compute the terms d\J\ and <9 2 J 2 separately. First, d\ J\ becomes 



dMt) =limtrT 3 r 1 



limtr T3Y1 



iM-n %M-U 

i i 



lim tr r 3 



s— >t 



iM-n iM-n 

1 1 



7iH<9D 



iM-n iM-n 



5 2 (t-s) 
5 2 (t - s) 



S 2 (t-s) 



(A-2) 



Next, for 82 J 2, we use eq. (15. 9p . 



d 2 J 2 (t) =-hmtrr 3 r 2 ^ 7 W [{-id\) n - — ^5 2 (t - s) 

iM -n 



n>l 



limtrr 3 r 2 y^7 : |. n) 



lim tr r 3 



n>l 



-idl) 



1 



72 {-idl) 



iM-n iM-n 
1 1 



[id s 2 r 



5 2 (t-s) 



iM-n iM-n 



l2{-idl) 



5 2 {t-s). (A-3) 
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Taking these into account, we obtain 



djJj(t) = limtrT3 



(sn - U) — + -^(6H - H) 

iM-n %M-n 



5 2 {t-s). (A-4) 



Note that {5H, r 3 } = 0. Thus, we obtain 



djJJt) = -21imtrr 3 =8\t - s). (A-5) 

Appendix B: Detailed calculation of J 2 

In the definition of J 2 in eq. (15. 6p . T 2 in eq. (15. 7p is not a simple matrix but includes 
complicated derivatives. We show in this Appendix that in the plane wave basis, it 
can be written simply as the derivative of H B with respect to k x , as in eq. (|5.27j) . 
Representing the 8 function as the plane wave, we see that 

J 2 (k v ,±oo) =limtrT3r 2 —S 2 (t — s) 

iM-Ti 

= l im f ^L t TT 3 T 2 = - -e m ^- s l (B-l) 

-*t] (2vr) 2 2 iM + in [vd^-nn^ky)] 1 ; 

Using (15.81) . we have 

T 2 e iU ^ = e in ( t -»)r 2 ^ 7 Wp n (fc s - id\, -k x - id s 2 ) 

n>l 

= e in ^r 2 J2li n) Pn(k x , -k x ) + 0{-id 2 ). (B-2) 

n>l 

Here, the derivative terms can be neglected, since n B does not depend on x. Now, note 
that 

Pn(k x ,-k x )= H^r'"^-*.)' = nkr 1 . (B-3) 

0<i<n-l 

Then, we obtain 

81-/ B 

E ^Pnfe' -*.) = E 1 = (b-4) 

n>l n>l x 

Therefore, after taking the M — > limit as well as s — > t, we have 

n f d2fi /0H B \ 1 

J (2tt) 2 V dk x J n [tvun + vd ky - r 3 n B \ 

Appendix C: Calculation of the Chern number 

We show that the r.h.s. of eq. fl5.3Uj) is simply the Chern number in eq. (I3.4p . Note 
that the former equation can be written as 

1 f dk y J 2 (k y ) = ^ / d 3 Mr G 2 (^Pj G, (CI) 



2 J y zv y > (2tt) 2 J \dk x J \dk 
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where d 3 k = dk y dQ = d 2 kdu. Inserting a partition of unity, 

X>«(*#i(*) = 1, (C2) 

Q 

as well as 

ipf}(d ki H B )ij a = (d kl £a)S/3a - (e/3 - e a )i)\{d ki i) a ) : (C-3) 

where e a (fc) and ip a (k) are respectively the eigenvalue and corresponding eigenfunction 
of the bulk Hamiltonian defined in eq. (13. II) . we obtain 

1 

2 



(2tt) 2 

The integration over a; yields 

^ I dk y J 2 {k y ) = ^- I d'k J2 [(dk^MUdkja) ~ (h k y ^ . (C'5) 

It turns out that this is simply the Chern number computed in eq. (13.41) . 
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